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We  have  used  tools  from  theory  of  harmonic  analysis  and  number  theory  to  extend  existing  theories  and  develop  new  approaches  to  problems.  This 
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work.  We  are  also  working  on  extending  our  work  to  multiply  periodic  processes.  We  have  also  extended  our  work  on  multichannel  deconvolution. 
We  have  used  the  tools  from  multichannel  deconvolution  to  develop  a  new  procedure  for  multi-rate  sampling.  We  have  investigated  applying  these 
techniques  to  develop  a  new  procedure  for  A-D  conversion.  We  have  also  developed  interlinked  wavelet  bases,  interlinked  via  number-theoretic 
conditions  that  proved  useful  for  both  multichannel  deconvolution  and  multi-rate  sampling.  We  have  also  created  new  procedures  for  sampling  in 
radial  domains,  exploiting  coprime  relationships  among  zeros  of  Bessel  functions.  These  results  have  application  to  radar  and  sonar.  We  have  also 
extended  these  ideas  to  operator  theory,  creating  sets  of  strongly  coprime  chirp  and  chirplet  operators.  Finally,  we  have  developed  a  new  way  of 
visualizing  mappings  of  the  complex  plane.  We  have  a  given  function  evolve  in  time,  starting  from  the  unaltered  complex  plane  and  ending  with  the 
range  of  the  function.  We  use  a  computer  to  develop  a  frame-by-frame  movie  of  the  evolution,  which  appears  to  be  a  continuous  deformation.  We 
are  using  this  technique  to  create  an  electronic  dictionary  of  conformal  mappings. 
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1. )  FORWARD 

We  have  used  tools  from  theory  of  harmonic  analysis  and  number  theory  to  extend  existing  theories 
and  develop  new  approaches  to  problems.  This  work  has  focused  on  several  areas.  We  have 
developed  algorithms  on  extensions  of  Euclidean  domains  which  have  led  to  new  computationally 
straightforward  algorithms  for  parameter  estimation  for  periodic  point  processes,  and  in  particular, 
for  sparse,  noisy  data.  This  is  useful  in  the  analysis  of  radar  and  sonar  data,  among  other  things.  We 
have  shown  why  Fourier  analytic  methods,  e.g.,  Wiener’s  periodogram,  do  not  produce  maximum 
likelihood  estimates  for  the  sparse  data  sets  on  which  our  methods  work.  We  are  also  working  on 
extending  our  work  to  multiply  periodic  processes.  We  have  also  extended  our  work  on  multichannel 
deconvolution.  We  have  used  the  tools  from  multichannel  deconvolution  to  develop  a  new  procedure 
for  multi-rate  sampling.  We  have  investigated  applying  these  techniques  to  develop  a  new  procedure 
for  A-D  conversion.  We  have  also  extended  the  work  on  both  deconvolution  and  sampling  to 
radial  domains,  exploiting  coprime  relationships  among  zero  sets  of  Bessel  functions.  We  have  also 
discussed  applications  of  these  ideas  to  specific  applied  problems,  including  radar  and  sonar.  We 
have  developed  interlinked  wavelet  bases,  interlinked  via  number-theoretic  conditions  that  proved 
useful  for  both  multi-channel  deconvolution  and  multi-rate  sampling.  We  have  extended  these 
ideas  to  operator  theory,  creating  sets  of  strongly  coprime  chirp  and  chirplet  operators.  Finally, 
we  have  developed  a  new  way  of  visualizing  mappings  of  the  complex  plane.  We  have  a  given 
function  evolve  in  time,  starting  from  the  unaltered  complex  plane  and  ending  with  the  range  of 
the  function.  We  use  a  computer  to  develop  a  frame-by-frame  movie  of  the  evolution,  which  appears 
to  be  a  continuous  deformation.  We  are  using  this  technique  to  create  an  electronic  dictionary  of 
conformal  mappings. 
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3.)  STATEMENTS  OF  PROBLEMS  STUDIED 


Our  work  on  multichannel  deconvolution  included  the  following. 

i. )  Fit  our  general  models  of  system  kernels  with  models  of  currently  deployed  active  and  passive 
remote  sensors  as  convolution  equations.  Then  apply  our  theory  to  develop  systems  in  which 
complete  signal  information  can  be  recovered.  In  particular,  we  are  looking  to  apply  the  theory  in 
general  imaging  systems  (e.g.,  for  sub-pixel  resolution),  radar  (e.g.,  strongly  coprime  chirp  pulses), 
etc.  (see  [2],  [8],  [38]). 

ii. )  Develop  a  basis  in  which  to  develop  the  multichannel  theory  which  is  more  amenable  to 
discretization.  This  would  then  be  directly  applicable  to  (i.)  (see  [2],  [8]). 

iii. )  Coordinate  the  theory  with  filtering  systems  for  the  removal  of  noise  and/or  other  unwanted 
information  (see  [15]). 

iv. )  Coordinate  the  theory  with  irregular  sampling  theory,  and  wavelet  and  Gabor  analysis  (See 

[11],  [12],  [10],  [2]). 


We  also  continued  our  work  on  parameter  estimation. 

v. )  Develop  computationally  straightforward  techniques  for  spectral  analysis  of  a  very  broad  class 
of  periodic  processes,  including  procedures  so  that  estimates  achieve  the  Cramer-Rao  bound  (see 

[9])- 

vi. )  Extend  these  techniques  to  the  complete  analysis  of  multiply  periodic  point  processes,  including 
the  recovery  of  the  fundamental  period(s),  phase  information,  the  multiples  of  the  periods,  and  the 
deinterleaving  of  the  data  (see  [7]). 

Additional  items  include  new  results  on  derivatives  and  nronotonicity  (see  [5] ,  [6] )  and  new  proce¬ 
dures  for  visualizing  complex  mappings  (see  [3] ) . 


4.)  SCIENTIFIC  PROGRESS  AND  ACCOMPLISHMENTS 


Sampling  on  Unions  of  Non-Commensurate  Lattices  via 
Complex  Interpolation  Theory  [11],  [12] 

Main  Results 

Solutions  to  the  analytic  Bezout  equation  associated  with  certain  multichannel  deconvolution  prob¬ 
lems  are  interpolation  problems  on  unions  of  non-comnrensurate  lattices.  These  solutions  provide 
insight  into  how  one  can  develop  general  sampling  schemes  on  such  sets.  We  give  specific  examples  of 
non-comnrensurate  lattices,  and  use  a  generalization  of  B.  Ya.  Levin’s  sine-type  functions  to  develop 
interpolating  formulae  on  these  sets.  Let  t  E  M,  and  let  a  be  an  irrational  that  is  poorly  approxi¬ 
mated  by  rationals,  e.g.,  y/2,  or  (1  +  v/5)/2.  Let  Hi(t)  =  Xr_11](i)  ,  ^{t)  =  Xr_aja](i)  model  the 


2 


impulse  response  of  the  channels  of  a  two-channel  system.  Then  m  (()  =  blll(^7r^)  ?  /j,2(Q  = 
Since  a  is  poorly  approximated  by  rationals,  {m}  is  strongly  coprime.  Now  let 


id  = 


for  k  £  N.  and  T  =  Ti  U  T2  . 


Note  that  the  information  contained  in  the  original  signal  is  reconstructed  by  creating  deconvolvers 
defined  initially  on  T  U  {0}.  Order  the  elements  of  T,  denoting  this  as  T  =  {A/,.}.  We  have  the 
following. 


Theorem  1  Let  a  be  an  irrational  that  is  poorly  approximated  by  rationals,  and  let  f  be  a  (1  +  a)- 
band-limited  function.  Then  f  is  uniquely  determined  by  {/(A*,)}  U  {/(0),  f\ 0)}  . 

Note,  for  a  (l  +  Q!)-band-limited  function,  the  Nyquist  rate  is  l/(2(l+o:)).  However,  our  individ¬ 
ual  sampling  rates  are  1/2  and  l/(2a).  Both  these  rates  are  below  Nyquist.  The  reconstruction  of 
/  from  this  lattice  is  achieved  by  using  complex  interpolation  theory.  These  techniques  go  back  to 
basic  Lagrange  interpolation,  and  were  developed  for  entire  functions  by  various  mathematicians. 
We  first  solve  the  problem  on  TU  {0}.  Let  G(z)  =  sin(27rz)  sin(27raz)  .  Then,  G(z)  is  an  entire 
function,  which  is  almost  periodic  on  M,  and  has  simple  zeros  on  T,  and  a  double  zero  at  {0}. 


Proposition  1  Let  A&  G  F  and  let 

Hkiz)  =  awr-  m  '  andso  HiM  =  s»- 

At  z  =  0,  we  have  to  construct  interpolating  functions  K\ ,  iL2  so  that 

K\ (0)  =  1  ,  K[(0)  =  0  ,  K2{ 0)  =  0  ,  K'2{ 0)  =  1  . 
Using  the  Taylor  expansion  of  G,  we  derive  the  following. 


Proposition  2 

K  (r,  G{z)  G(z)  G(z )  G(z) 

11  j  (G"(0)/2!)(22)  (47 r2a)(z2)  ’  2[^  (G"(0)/2!)(z)  (4tt 2a)(z)  ' 

Combining  these  two  propositions  gives  us  the  reconstruction  formula. 


Theorem  2  Let  a  be  an  irrational  that  is  poorly  approximated  by  rationals,  and  let  f  be  a  (1  +  a)- 
band-limited  function.  Let  Ti  =  {^}  ,  T2  =  for  k  G  N,  and  Ietr1ur2  =  r  =  { A^}  . 

Then  f  is  uniquely  determined  by  {/(A*,)}  U  {/(0),  /'( 0)}  .  Moreover,  f  can  be  approximated  from 
its  values  on  T  U  {0}  by  the  formula 


m 


G(t) 


G'(\k)(t-\k) 


+  /( 0) 


G(t) 


(4:ir2a)(t2) 


+  a  o) 


G(t) 


(47 r2a)(t) 


where 

G(t )  =  sin(27rt)  •  sin(2a7rf)  . 
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We  need  to  point  out  two  items.  Let  r i  =  1,  ?’2  =  a.  First,  the  sampling  grid  is  rigid. 
Perturbation  of  the  grid  results  in  a  loss  of  information.  Second,  that  because  sampling  points  in  T 
can  get  arbitrarily  close  together  (inf{|Am  —  An|}  =  0,  and  so  inf{J][fc_^  |  sin(27r?’fcA)|  :  A  e  Tj}  =  0), 
the  dual  elements  of  the  interpolating  functions  can  not  form  a  Riesz  basis  and  the  interpolating 
formula  can  not  converge  in  norm.  In  fact,  the  dual  elements  of  the  interpolating  functions  do  form 
a  Bessel  sequence,  but  do  not  form  a  frame  and  therefore  do  not  form  a  Riesz  basis.  The  problems 
occur  at  points  where  the  sampling  points  get  close  together.  The  interpolating  function  follows  the 
original  function  along  exactly,  except  for  a  very  subtle  “ripple”  at  those  points  where  the  sampling 
points  get  close  together.  We  are  currently  exploring  ways  to  stabilize  the  construction.  The  first 
involves  continuously  backing  off  the  bandwidth.  This  works,  but  we  need  to  set  up  exact  bounds. 
The  second  approach  involves  rewriting  T  in  terms  of  the  “separated”  and  “close”  points.  At  the 
close  points,  we  may  be  able  to  rewrite  the  formulae  in  terms  of  f. 

The  result  generalizes.  We  can  create  sampling  sets  on  t  lattices  using  a  set  of  numbers  {rj}f=1 
such  that  (rt/rj)  is  poorly  approximated  by  rationals  for  i  /  j.  Again,  if  {pij^zl  is  a  set  of  primes, 

{!>  Vpi,  VP1P2,  ■■■,  VPW2  ■■■Pi- 1} 

is  a  set  of  numbers  whose  ratios  are  poorly  approximated  by  rationals.  Let  F/,. 
and  let 

l 

u  Ffc  =  r  =  {Afc} . 

k= 1 

We  reconstruct  on  V  U  {0},  letting 

tt  G(z) 

G(z)  =  TT  sin(27rrfcz)  and  letting  Hm(z )  =  *  . — - 

^=1  G\z)(z-  Am, 

Then  Hm( An)  =  5mn  . 

The  interpolating  functions  at  the  origin  are  a  linear  combination  of  gj  =  ,  j  =  1 , . . . ,  l , 

chosen  so  that  ^(0)  =  Skj  ,  k,j  =  We  can  also  create  these  as  unions  of  regular 

lattices  in  higher  dimensions. 

We  close  the  paper  by  proving  that,  from  a  measure-theoretic  and  thus  probabilistic  point  of 
view,  the  strongly  coprime  condition  on  which  our  constructions  lie  is  quite  natural.  The  numbers 
in  the  complement  of  our  condition  are  called  the  Liouville  numbers  L.  We  show  the  following. 

Theorem  3  The  class  of  Liouville  numbers  L  has  Lebesgue  measure  zero. 

Thus  if  one  chooses  an  irrational  number  at  random  (say,  using  the  uniform  distribution  on 
any  finite  interval) ,  one  would  almost  surely  choose  an  irrational  which  is  poorly  approximated  by 
rationals.  However,  the  situation  is  even  better  than  this. 

Theorem  4  Let  L  denote  the  class  of  Liouville  numbers.  Then,  for  every  a  >  0,  L  has  o- 
dimensional  Hausdorff  measure  zero. 


=  for  n  G  N, 


New  Directions  in  Sampling  and  Multi-Rate  A-D  Conversion 
Via  Number  Theoretic  Methods  [10] 
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Main  Results 


This  paper  uses  the  new  sampling  theory  developed  in  the  paper  above  to  provide  a  new  approach 
to  multi-rate  A-D  conversion.  This  uses  multiple  converters  with  different  rates.  Consider  a  1  +  a 
bandlimited  signal,  with  Nyquist  sampling  rate  fjy  =  2(1  +  a)  (samples  /  sec).  Let  the  sampling 
rates  be  denoted  fi,  i  =  1,  •  •  •  ,  K.  These  rates  are  such  that  (i)  their  ratios  are  strongly  coprime, 
and  (ii)  jr  fi  =  In,  (the  average  of  the  rates  is  equal  to  the  Nyquist  rate). 

As  above,  the  signal  may  be  reconstructed  from  two  samplers  running  at  rates  f\  =  2,  and 
/2  =  2a.  We  emphasize  that  this  choice  of  /i  is  arbitrary  and  may  take  any  value  so  long  as  the 
above  conditions  are  satisfied.  While  /i  remains  fixed  at  a  predetermined  rate,  /2  — ►  /v  from  below 
as  a  becomes  large.  Thus,  the  introduction  of  the  fixed  low  rate  sampler  at  f\  has  relaxed  the  faster 
sampling  rate  below  the  Nyquist  rate.  Because  the  multiple  rates  are  non-commensurate,  they  will 
eventually  yield  sampling  points  that  are  arbitrarily  close  to  one  another,  within  5  >  0  say.  In  this 
case,  there  exists  an  e  >  0  such  that  one  of  the  samples  occurring  within  the  5  neighborhood  may 
be  omitted,  while  guaranteeing  unique  reconstruction  for  signals  bandlimited  to  1  +  a  —  e.  It  is  also 
possible  to  look  at  reconstruction  f  locally. 

The  multi-rate  sampling  approach  can  be  extended  to  more  than  two  rates.  Consider  a  1  +  a+/3 
bandlimited  signal,  with  a  and  (3  selected  to  satisfy  the  coprime  condition  (for  example,  choose  a 
and  to  be  square  roots  of  products  of  primes).  Now,  three  parallel  samplers  may  be  used  with 
rates  fi  =  1/2,  /2  =  l/(2a),  and  fs  =  1/(2 (3).  (Again  we  emphasize  that  these  choices  are  not 
unique.)  We  consider  two  cases.  First,  if  a  «  fi,  then  /2  and  approach  /at/2  from  below  as  a 
increases.  Second,  if  /?  «  ka,  for  k  >  0  integer,  then  /2  — s ►  /jv/(1  +  A;),  while  fe  — >  /jyfc/(  1  +  k )  as 
a  grows.  Thus,  many  combinations  of  slow  and  fast  rates  are  possible. 

It  is  interesting  to  contrast  this  scheme  with  the  more  conventional  case  of  K  interleaved  iden¬ 
tical  parallel  A-Ds,  sampling  at  the  same  rate,  fl  =  /n/K-  The  coprime  multi-rate  approach  allows 
the  A-Ds  to  run  at  different  rates,  without  the  need  for  signal  preconditioning,  e.g.,  bandpass  fil¬ 
tering.  Thus,  the  proposed  multi-rate  scheme  generalizes  the  conventional  interleaved  case.  Also, 
the  slower  A-Ds  could  utilize  more  quantization  levels.  The  multi-rate  approach  requires  non- 
commensurate  rates,  so  a  scheme  must  be  developed  for  generating  these  rates  while  maintaining 
the  underlying  synchronism.  This  requires  close  approximation  to  irrationals  from  rational  num¬ 
bers.  In  addition,  the  multiple  rates  will  cause  sampling  times  to  become  arbitrarily  close.  However, 
as  noted  above,  when  neighboring  samples  fall  into  a  5  neighborhood,  one  sample  may  be  omitted 
(with  a  corresponding  e  small  reduction  in  the  achievable  signal  bandwidth). 


Sinusoidal  Frequency  Estimation  via  Sparse  Zero  Crossings  [9] 

Main  Results 

We  consider  estimation  of  the  period  of  a  sinusoid  in  additive  Gaussian  noise,  based  on  observations 
of  the  zero-crossing  (ZC)  times.  The  problem  is  treated  in  a  continuous-time  framework.  It  is 
assumed  that  the  signal-to-noise  ratio  is  sufficient  (roughly  >  8  dB)  such  that  the  noise  may 
be  approximated  as  additive  in  the  phase.  An  exact  mean-square  error  analysis  is  provided  for 
this  approximation.  We  apply  modified  Euclidean  algorithms  (MEA’s)  and  their  least-squares 
refinements  in  this  framework,  to  estimate  the  period  of  the  sinusoid,  with  low  complexity.  Unlike 
linear  regression  methods  based  on  phase  samples,  the  proposed  approach  works  with  very  sparse 
ZC  measurements,  and  is  resistant  to  outliers.  The  MEA-based  approach  is  motivated  by  the  fact 
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that,  in  the  noise- free  case,  the  greatest  common  divisor  (gcd)  of  a  sparse  set  of  the  first  differences 
of  the  zero  crossing  times  is  very  highly  likely  to  be  the  half-period  of  the  sinusoid.  The  MEA  acts 
to  robustly  estimate  the  gcd  of  the  observed  noisy  data.  The  MEA  period  estimate  may  be  refined 
via  a  least-squares  approach,  that  asymptotically  achieves  the  appropriate  Cramer-Rao  bound. 
Simulation  results  illustrate  the  algorithms  with  as  few  as  10  zero-crossing  times.  The  algorithm 
behavior  is  also  studied  using  Bernoulli  and  random  burst  models  for  the  missing  ZC  times,  and 
good  performance  is  demonstrated  with  very  sparse  observations. 

The  observed  data  is  modeled  as  a  time  series  -  a  sampled  sinusoid  in  Gaussian  noise,  expressed 
as  x(n)  =  A  exp  i(u>onp  +  f>)  +  z(n)  ,  for  n  =  no,  no  +  1, . . . ,  no  +  N  —  1.  The  noise  z(n)  is  assumed 
to  be  a  white  Gaussian  zero-mean  complex  noise  sequence  with  variance  a2.  The  amplitude  A, 
frequency  u>q,  and  phase  are  regarded  as  unknown  constants  for  a  particular  realization.  The 
value  p  is  the  sampling  period.  We  can  write 


x{n)  =  Aexpi(w>onp  +  <f)  +  z(n)  = 

=  [1  +  n(n)]  Aexp(i((u0np  +  0))  , 


Z\  77, ) 

1  H - r -  exp(-i(LUonp  +  <f>)) 


A 


Aexp(i(u;onp  +  4>)) 


where  v(n)  =  exp  [— i(w>onp  +  4>)\.  Note  that  E[i/(n)\  =  0  and  var{n(n)}  =  cr2/A2  =  SNR'1,  so 

u{n)  — >  0  as  SNR  — »  oo.  If  we  split  v  into  real  and  imaginary  parts,  writing  v  =  w  +  iu,  we  can 
express  1  +  n(n)  in  terms  of  modulus  and  phase  as 


1  +  u(n) 


[1  +  w(n )]2  +  [n(n)]2 


exp 


arctan 


u(n)  \ 
1  +  w(n)  ) 


We  give  two  proofs  of  the  following  result  in  the  paper. 


Theorem  5  If  A2 /a2  2>  1,  (i.e.,  high  SNR),  then  [1  +  v{n)]  ~  exp  [iu(n)\  .  This  in  turn  gives 
x(n)  =  Aexpi(uonp  +  <f)  +  z(n)  ~  Aexp(i(uionp  +  0  +  u(n)))  . 


Residue  and  Sampling  Techniques  in  Deconvolution  [8] 

Main  Results 

In  this  paper  we  will  study  techniques  for  solving  the  following  general  problem,  which  we  will  refer 
to  as  the  Multisensor  Deconvolution  Problem  or  MDP:  Given  a  collection  of  compactly  supported 
distributions,  1  C  how  can  we  recover  an  arbitrary  function  f  6  C'°°(Md)  from  the 

data  {sj}^  =  {/  *  The  solution  technique  that  we  consider  involves  the  construction  of 

deconvolvers,  which  come  in  a  variety  of  types  but  which  are  essentially  a  collection  of  distributions 
which  (1)  depend  only  on  the  convolvers  {/u}™  i  and  (2)  allow  for  the  solution  of  the  MDP  with 
only  simple  linear  operations  on  the  data  {s,}^. 

We  describe  a  solution  to  a  modified  analytic  Bezout  equation.  The  technique  of  solution 
involves  the  Jacobi  interpolation  formula  and  the  Cauchy  residue  calculus.  We  create  an  approxi¬ 
mate  deconvolution  system  by  solving  the  Bezout  equation  Xa=i  hi '  =  Ti  where,  for  0  <  e  <  1, 

<p(x)  =  ipe{x)  =  e  1  'f>(e~1x).  Thus,  — >  1  as  e  — >  0+.  This  gives  us  our  deconvolving  functions. 
Then,  as  ip  — >  5,  f  *  <p  — >  /  in  the  sense  of  distributions.  That  is,  the  signals  st  are  filtered  by 
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the  vl  tr,  (which  have  been  created  digitally,  optically,  etc.,  in  coordination  with  the  creation  of  the 
system  and  possibly  tailored  to  be  optimized  under  some  constraint)  and  added,  resulting  in  the 
reconstruction  of  /  *  ip.  By  controlling  e,  /  can  be  reconstructed  to  any  predetermined  accuracy. 
We  also  create  these  solutions  by  using  techniques  from  real-variable  sampling  theory. 


On  Positive  Derivatives  and  Monotonicity  [5],  [6] 

Main  Results 

A  fundamental  result  of  calculus  states  that  a  function  /  with  a  positive  derivative  Df  on  an 
interval  is  increasing  on  that  interval.  The  result  follows  directly  from  the  Mean  Value  Theorem. 
We  explore  the  extent  to  which  the  hypotheses  of  the  Mean  Value  Theorem  can  be  weakened  and 
/  still  shown  to  be  increasing.  Let  /  be  a  continuous  function  on  an  interval.  By  constructing 
counterexamples  using  Cantor  sets,  we  show  that  the  assumption  Df  >  0  a.e.  does  not  imply 
that  /  is  increasing.  We  can,  however,  show  that  if  Df  >  0  except  on  a  countable  subset  of  an 
interval,  then  /  is  increasing.  We  call  this  the  Countable  Exceptional  Set  Theorem.  This  theorem 
is  generalized  by  the  Goldowsky-Tonelli  Theorem,  which  tells  us  that  if  Df  exists  except  on  a 
countable  subset  of  an  interval  and  Df  >  0  a.e.,  then  /  is  increasing.  However,  we  then  go  on 
to  show,  in  a  very  natural  sense,  that  Goldowsky-Tonelli  is  a  vacuous  extension  of  the  Countable 
Exceptional  Set  Theorem.  We  close  by  exploring  a  few  related  questions,  asking  to  what  extent 
nronotonicity  implies  positive  derivative. 


Complex  Analysis  An  Animated  Approach  [3] 

Sample  Lessons 

Sample  Lesson  1:  An  Introduction  to  Animated  Graphs 

One  of  the  most  fundamental  skills  needed  in  understanding  complex  analysis  is  visualizing  complex 
functions.  This  is,  however,  a  challenging  task.  The  graph  of  a  complex  function  /  :  C  — >  C  lives 
in  CxC  (two  complex  dimensions),  which  is  four  real  dimensions.  Given  that  we  only  live  in  three 
spatial  dimensions,  visualization  is  difficult.  We  can,  however,  get  around  this  by  using  our  fourth 
dimension  -  time  -  as  a  tool  to  help  us  visualize  complex  functions.  We  will  watch  a  given  function 
evolve  in  time,  starting  from  the  unaltered  complex  plane  and  ending  with  the  range  of  the  function. 
The  computer  provides  us  with  an  ideal  tool  to  see  this  evolution.  We  compute  a  frame-by-frame 
movie  of  the  evolution  “off-line,”  and  show  it  at  relatively  high  speeds,  giving  the  appearance  of  a 
continuous  deformation.  These  “movies”  of  complex  functions  will  be  called  animated  graphs. 

A  mathematical  framework  for  animated  graphs  already  exists.  This  deformation  of  the  function 
is  a  homotopy,  from  the  identity  to  the  target  function.  We  will  give  a  formal  definition  of  homotopy 
elsewhere.  For  now,  we  can  think  of  these  graphs  as  movies  in  time  t,  0  <  t  <  T,  where  we  are 
looking  at  F(z,  t)  =  (1  —  t)z  +  tf(z)  .  Since  a  complex  function  could  have  a  domain  as  large  as  the 
entire  complex  plane,  we  will  highlight  some  of  the  interesting  curves  and  arcs  in  the  domain. 

The  animated  graphs  were  developed  using  f(z),  a  software  package  for  complex  variables  written 
by  Martin  Lapidus  of  Lascaux  Graphics.  The  animated  graphs  are  programs  developed  by  the 
author  within  f(z),  just  as  one  would  write  Notebooks  inside  of  Mathematica.  The  animated  graphs 
consist  of  curves  and  arcs  in  the  domain  on  the  left  and  their  images  under  mapping  on  the  right. 
However,  it  is  very  important  to  note  that  you  will  not  just  see  the  static  images,  but  rather  a 
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“continuous”  (actually,  frame-by-frame)  deformation  of  the  curves  and  arcs  in  the  domain  to  their 
corresponding  images  in  the  range.  It  is  also  important  to  note  that  this  is  an  extension  of  the 
traditional  representation  of  complex  mappings,  namely  by  drawing  curves,  arcs,  and  or  regions  in 
the  domain,  and  then  showing  their  images  in  the  range. 

We  also  have  a  set  of  animated  graphs  for  linear  mappings,  and  will  expand  upon  the  matrix 
model  of  C.  The  animated  graphs  will  include  linear  mappings,  and  will  demonstrate  that  the 
action  of  multiplying  by  a  complex  number  is  the  same  as  the  action  of  a  2  x  2  matrix  on  the 
vector  z.  The  animations  of  the  linear  maps  show  the  actions  of  the  maps  on  rectangular  grids. 
For  example,  the  mapping  iz  =  em/2rel6  =  re^+^/ 2  ro^a^es  a  COnrplex  number  counterclockwise  by 
7t/2.  Note  that  this  is  the  same  result  as  multiplying  the  complex  vector  z  =  (x,  y )  by  the  matrix 


Mt  = 


-1 

0 


Sample  Lesson  2:  Local  Behavior  of  Polynomial  and  Rational  Functions 

This  second  group  of  sample  lessons  is  included  to  show  that  the  animated  graphs  can  be  used 
to  effectively  demonstrate  certain  principles.  An  important  concept  in  complex  variables  is  that, 
thanks  to  the  power  of  Taylor  and  Laurent  expansions,  the  behavior  of  an  analytic  or  meromorphic 
function  in  a  neighborhood  of  a  zero  or  a  pole  can  be  understood  by  factoring  the  function.  For 
example,  let  /  be  an  analytic  function  which  has  a  zero  of  order  k  at  zq.  In  a  sufficiently  small 
neighborhood  A f  of  zq,  f  can  be  written  as  f(z)  =  (z  —  zo)kF(z)  ,  where  F  is  analytic  and  non-zero 
in  N .  Moreover,  the  behavior  of  the  mapping  /  in  AA  is  essentially  f(z)  «  (z  —  zo)kF(zo) ,  that 
is,  the  mapping  looks  like  the  polynomial  (z  —  zo)k  multiplied  by  the  complex  number  F(zq).  This 
important  principle  is  demonstrated  for  the  functions  f(z)  =  (z  —  (l  +  i))(z  —  (^  —  and  f(z)  = 

z4  —  1.  The  animations  show  how  the  function  maps  neighborhoods  of  the  zeros  onto  a  region 
in  a  neighborhood  of  the  origin  in  the  range.  Even  though  the  circles  and  line  segments  in  the 
neighborhoods  of  zeros  end  up  on  the  same  curves  in  the  image,  the  paths  they  take  to  get  there 
are  quite  different.  The  different  magnifications  and  rotations  that  the  circles  and  line  segments 
undergo  are  the  result  of  multiplying  the  map  (z  —  zq )  by  the  appropriate  F(zq). 

Similarly,  let  g  be  a  meromorphic  function  which  has  a  pole  of  order  l  at  wq.  In  a  sufficiently 
small  neighborhood  Ad  of  wg ,  g  can  be  written  as  g{z)  =  (z  —  wg )~lG(z)  ,  where  G  is  analytic  and 
non-zero  in  Ad.  Moreover,  the  behavior  of  the  mapping  5  in  Ad  is  essentially  g(z)  ~  (z— wg)~lG(wg)  , 
that  is,  the  mapping  looks  like  the  function  (z  —  wg)~l  multiplied  by  the  complex  number  G(wg). 
This  important  principle  is  demonstrated  for  the  functions  g(z)  =  ,_^+^  and  g(z)  = 

Let  /  be  a  non-trivial  analytic  function  in  a  region  11  with  a  zero  of  order  k  at  zg  £  If.  Thus  by 
the  Taylor  expansion  of  /,  f(z)  =  Y,n=k  n!^  (z  ~  zo)n  ,  because  f(z0)  =  f(zQ )  =  f"(zQ)  =  ...= 
/(fc_1)(2o)  =  0  .  Also,  since  zeros  of  a  non-trivial  analytic  function  may  not  cluster,  there  exists  a 
neighborhood  M  in  which  we  may  write  f(z)  =  ((z  —  zo)k)F(z)  ,  where  F  is  a  non-zero  analytic 

function  in  M  such  that  lhnz^Zo  F(z)  =  -  /  0  .  We  get  the  Taylor  series  for  F  by  factoring 

F(z)  =  n\Z° ^  (z  ~  zo){n~k)  ■  Therefore,  in  Af,  the  behavior  of  the  mapping  /  is  essentially 

f(z)  ~  (z  —  zo)kF(zo) ,  that  is,  the  mapping  looks  like  the  polynomial  (z  —  zo)k  multiplied  by  the 
complex  number  F(zq). 

We  will  first  demonstrate  this  important  principle  for  the  function  f(z)  =  (z— (l+i))(z— (^— i^)) 
has  two  simple  zeros.  In  a  neighborhood  of  1  +  i,  f(z)  ~  +  i|)(z  —  (1  +  i))  =  f'(l  +  i)(z  —  (1  +  *))  • 

In  a  neighborhood  of  (|  —  i\),  f(z)  «  +  iFf){z  —  —  i|))  =  f'{\  —  i\){z  —  (|  —  i\))  . 


The  function  f(z)  =  z4  —  1  has  a  little  more  interesting  behavior.  Neighborhoods  of  each  of  the 
four  simple  zeros  l,i,  —1,  —i  get  mapped  to  neighborhoods  of  the  origin.  However,  each  of  these 
neighborhoods  takes  a  different  path  to  get  there.  We  can  compute  the  different  rotations  and 
magnifications  that  each  neighborhood  undergoes  by  factoring.  If  we  index  the  zeros  as  z\  =  1, 
^2  =  i,  z3  =  -1,  24  =  -i,  then  f(z)  =  (z-  Zi)F(zi)  =  (z-  Zi)f(zi)  ,  where  F(z1)  =  (21)3  +  (21)2  + 
(21)  +  1  =  4  =  f'(zi) ,  F(z2)  =  (z2)3  +  i{z2)2  -  (22)  -  i  =  ~4 i  =  f(z2 ) ,  F(z3)  =  (23)3  -  (z3)2  + 
(23)  —  1  =  — 4  =  f(z3) ,  +(24)  =  (24)3  —  i{z^)2  —  (24)  +  i  =  4 i  =  /'(z^)  .  This  same  behavior  extends 
to  general  analytic  functions. 

Similarly,  let  g  be  a  meromorphic  function  which  has  a  pole  of  order  l  at  wq.  Using  the 

Laurent  expansion  in  a  deleted  neighborhood  JA  of  wq,  we  have  g(z)  =  an(z  ~  wo)n  + 

El=i  bm(z  -  ,  where  the  coefficients  are  given  by  the  integrals  an  =  2+I  Jr  (z-w*)n+i  dz , 

bm  =  2^  fr  ff(z)(z—wo)m~1  dz  ,  for  T  =  {2  :  \z—wq\  =  r}  contained  in  M.  Of  course,  we  can  use  the 
techniques  given  in  section  6.4  to  compute  the  coefficients.  We  can  once  again  factor  this  expansion, 
writing  g  as  g(z)  =  ((2  —  wq)^1^)G(z)  ,  where  G  is  a  non-zero  analytic  function  in  M  such  that 
lim~^0  G(z)  =  bk  ^  0  .  The  Taylor  series  for  G  is  G(z)  =  Y!m=i  b™(z  ~  +  J2™=o  an{z  ~ 

wo)(n+0  •  Therefore,  in  A4.  the  behavior  of  the  mapping  g  is  essentially  g(z)  &  (2  —  wq )  1G(wq ) , 

that  is,  the  mapping  looks  like  the  rational  function  (2  —  multiplied  by  the  complex  number 

G(w0). 

This  important  principle  is  demonstrated  for  the  function  g(z)  =  zrzi  >  which  has  four  simple 
poles  at  i,i,  —1,  —i.  We  can  compute  the  different  rotations  and  magnifications  that  each  neigh¬ 
borhood  undergoes  by  factoring.  If  we  index  the  poles  as  w\  =  1,  w2  =  i,  w3  =  —1,  uq  =  —i,  then 
9(z)  =  ^±^G(Wi) ,  where  G(w i)  =  {wi)3+{w^+{wi)+1  =  \  ,  G(w2)  =  (^;2)3+,(^21)2-(^;2)-^  =  ^  = 

4  ’  G(w%)  (1^3)^  —  (iy3)a+(iy3)— 1  4  ’  ^(^4)  {w^)‘6  —i(w^y2  —  (1^4) 4 i  4 

This  same  behavior  extends  to  general  meromorphic  functions.  However,  the  situation  is  quite 
different  in  a  neighborhood  of  an  essential  singularity.  Recall  that  the  Casorati-Weierstrass  Theorem 
told  us  that  a  deleted  neighborhood  of  an  essential  singularity  gets  mapped  densely  onto  the  entire 
complex  plane.  The  behavior  is  actually  more  amazing.  One  of  the  theorems  we  will  prove  is  the 
Great  Picard  Theorem,  which  tells  us  that  if  a  function  g  is  analytic  in  a  deleted  neighborhood 
of  a  point  wq  and  has  an  essential  singularity  at  wq  ,  then  in  every  deleted  neighborhood  of  iuq,  g 
assumes  every  complex  number,  with  one  possible  exception,  an  infinite  number  of  times. 

Systems  of  Convolution  Equations,  Deconvolution,  Shannon  Sampling, 
and  the  Wavelet  and  Gabor  Transforms  [2] 

See  enclosed  table  of  contents. 
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